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Abstract 

We present results on co-recursive associated Laguerre and Jacobi polynomials 
^ ' which are of interest for the solution of the Chapman-Kolmogorov equations of some 

. birth and death processes with or without absorption. Explicit forms, generating func- 

tions, and absolutely continuous part of the spectral measures are given. We derive 
t — ' fourth-order differential equations satisfied by the polynomials with a special attention 

■ to some simple limiting cases. 
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: 1 Introduction 



Starting from a sequence of orthogonal polynomials {Pn}n>o defined by the recurrence rela- 
tion 

Pn+2{x) = {x - Pn+i)Pn+i{x) - 7„+iP„(a;), n>0, (1) 
and the initial conditions 

Po(x) = l, Piix)=x-Po, (2) 
with 7„ G C and 7„ 7^ 0, several modifications were considered: 

• Associated polynomials arise when we replace n by n + c in the coefficients /3„ and 7^ 
(keeping 7„ 7^ 0). If c is an integer k these polynomials are called associated of order 
k. The associated polynomials of order one are the numerator polynomials. 

• Co-recursive polynomials arise when we replace /3o by /3o + /i. 
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• Perturbed polynomials arise when we replace 71 by A71, (A > 0). 

• Co-recursive of this perturbed polynomials arise when the two previous modifications 
are made together. 

• Generalized co-recursive and perturbed polynomials arise when we change /3„ and/or 
7„ at any level n. 

In the study of birth and death processes orthogonal polynomials, in particular all the 
hypergeometric families of the Askey scheme |]1|, and their corresponding associated fam- 
ilies, play a primordial role in the Karlin-McGregor solution of the Chapman-Kolmogorov 
equation [1^, |TT| 

Pmn{t) ~ / e-"*P^(x)P„(x)rf0(x). (3) 







In certain birth and death processes, zero-related polynomials []T2|, arise in a natural 

way. Zero-related polynomials are special co-recursive associated polynomials. 

More generally co-recursive and generalized co-recursive polynomials are involved in the 
solution of the Chapman-Kolmogorov equation of birth and death processes with absorption 
or killing 0,|T7|,|T§. 

The purpose of this paper is to present some results on co-recursive associated Laguerre 
and Jacobi polynomials which are of special interest in the study of birth and death processes 
with linear and rational rates respectively. The Laguerre polynomial families are involved in 
processes for which the birth and death rates are of the form 



Xn = n + a + c+l, fXn+i = n + c + 1, n > 0, /io = c - /i. (4) 

The Jacobi polynomial families are involved when the rates are of the form 

^ _ 2{n + c + a + i3 + l){n + c + p+l) ^ ^ ^ 

{2n + 2c + a + l3 + l){2n + 2c + a + f3 + 2y ' 

2{n + c) {n + c + a) . . 

" {2n + 2c + a + (3){2n + 2c + a + (3 + l)' ^ ' ^ ^ 

_ 2c(c + a) . 

~ {2c + a + p){2c + a + (3 + l) ~ ^' 



In both cases fio = correspond to the "honest" ||2^ linear processes (i.e. processes for 
which the sum of the probabilities Pmnif) is equal to 1). Cases /io = Const. 7^ correspond 
to processes with absorption and are not "honest" . However if /iq = c in the Laguerre case 
or /io = i^2c+a+p){2c+a+p+i) Jacobi case, the corresponding processes are simply solved 

using associated polynomials. 

In section ^ we explain the method used by applying it to the Laguerre case. In we 
give an explicit expression for the co-recursive associated Laguerre (CAL) polynomials. In 
|2?^ we derive a generating function of them and we found the absolutely continuous part 
of the spectral measure. The subsection |2]^ is devoted to the derivation, using the Orr's 



method, of a fourth-order differential equation satisfied by the CAL polynomials. In O we 
present results obtain in some limiting cases among which a new simple case of associated 
Laguerre polynomials. 
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In section ^ we present briefly some results corresponding to the Jacobi case. In ^Tl] we 
give an explicit expression for the co-recursive associated Jacobi (CAJ) polynomials. In |3.2| 
we present a generating function and in p.3| we give the absolutely continuous component 



of the spectral measure. The subsection is devoted to some limiting cases of CAJ poly- 
nomials for which we give fourth-order differential equations they satisfy. In section p.6| are 
some concluding remarks. 

We use the notation of for the special functions used in this work. We don't give 
validity conditions on the parameters of the used hypergeometric functions, analytic contin- 



uations or limiting processes giving, in general, valid formulas. We use Slater's notation p9 
for the product of T functions 

2 The case of Laguerre polynomials 

Replacing n by n -|- c in the recurrence relation of the Laguerre polynomials we obtain the 
recurrence relation satisfied by the associated Laguerre polynomials L'^{x] c) 

{2n + 2c + a + 1 - x)p„ = (n + c + l)p„,+i + {n + c + a)p„_i. (9) 

To complete the definition of the polynomials L'^{x] c) the initial condition 

L-_,{x;c) = 0, L-(x;c) = l, (10) 

has to be imposed. These polynomials are orthogonal with respect to a positive measure 
when {n + c){n + a + c) > 0, Vn > 0. (See for details). 

Note that if we consider the monic polynomials L"(x; c) defined by 

L^(x;c) = (-l)"(c+l)„L^(x;c) (11) 

they satisfy the recurrence relation 

{x-2n-2c-a- l)L°(x; c) = L^+i(a;; c) + {n + c){n + c + a)L'^_^{x; c). (12) 

We can see that this recurrence is invariant in the transformation T defined by 

T(c, a) = (c + a, — a). (13) 

2.1 Explicit representation of the CAL polynomials 

Equations (||) and (|l^) give for L"{x; c) 

L"(x;c) = —(x-2c-a-l) (14) 

c + 1 

The CAL polynomials L'^{x]c,fi) satisfy the same recurrence relation (|^) with a shift fi on 
the monic polynomial of first degree, i.e. 

L1{x]c,fi) = —^^{x + fi-2c-a-l). (15) 
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To obtain (x; c, /x) with (|^) we have to impose the initial condition for the CAL polynomials 



L'^iix; c, fi) 



fx 



c + a 



Lq{x; c,fi) = 0. 



(16) 



Even for c + a — this initial condition in the recurrence relations (^) leads to (|T5D . 

We know two linearly independent solutions of (§) in terms of confluent hypergeometric 
functions: 



{c + a + l)n „ -n - c 
(c+l)n \l + a 



and 



-n — c — a 
1 — a 



(17) 



Writing the polynomials L'^{x; c, fi) as a linear combination 

L"{x; c, fi) = Aun + Bvn 
and using the initial condition (|l^) we obtain 



A 



A 



c + a 



vo - v^i 



and 



B 



1 

'a 



c + a 



Uo - U-i 



(19) 



where A may be calculated using contiguous relations of confluent hypergeometric functions 
[|, page 253-254] 



A = M_ifo — UqV^i 



a 



c + a 



With the help of the relation 
^1-7 



7iFi 



X 



a 



(3^FA =(7-/3)2i^2 



'-7,/9-7 + l 
, a,/?-7 



(20) 



(21) 



we can write the CAL polynomials as 



L"(x; c, fi) 



(C+ l)n 

Ji — c 



l + T) 



X ■ 



, . /-c, /i - c + 1 \ /-n - c - a 



(22) 



This representation is the one we use to derive the fourth-order differential equation in 
section p.3| . It is valid only for a 7^ 0, ±1, ±2, . . . but these restrictions can be removed by 
limiting processes. 

Following the same way as in we can find an explicit representation. We first transform 
the iFi in (|r^) using Kummer's transformation page 253] 



X 



'c — a 



-X 



(23) 



Then we use formula 



iFiL;xUFr 



—X 



E 

fc=0 



p -k, l-k-b,c 

k\{h)k \ l-k-a,d '\ 



(24) 
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for the four products in (|T^). 

Applying the three term relation Eq. (2) page 15] to the four resulting 3^2(1) gives a 
sum of two terms which we can group to obtain the explicit form 



n\ ^ (c+ l)fc(c + a + l)fc 

x,fJ ^-^:"'" + "'^ + ? (25) 

^■^\c + k + l,c + a + k + l,G'J ' ^ ^ 

where 

G = (26) 

^^ 

Take care of the limiting processes when n = k and when c or c + a = 0. The representation 
can also be proved using the generating function (RO). 



2.2 Generating function, spectral measure 

Let F{x,w) be a generating function of the CAL polynomials L^{x; c, fi) 

00 

F(x,t/;) = ^^/;"L^(x;c,/i). 

n=0 



(27) 



The recurrence relation (^) and the initial condition (|T6D lead to the following differential 
equation for F{x, w) 

d 

w{l — wY-^F{x, w) + [(1 — w)(c — (c + a + l)w) + xw]F{x, w) = c — fxw. (28) 

The function F{x, w) is normalized by the condition F{x, 0) = 1 and due to the orthogonality 
of the L^{x; c, /i), we have the boundary condition 



F{x, w)d(f){x) = 1, 



(29) 



where d(j){x) is the spectral measure. 

The solution of the differential equation ( pSD which is bounded at = is easily obtain, 
for c > 0, following the same methode as in ||12[ 



F(x,w) = w ''(1 — w) 



-a-l 



exp 



X 



X 



1 — w 
Jo 



1-u 



c — fiu) exp 



X 



l-u 



Changing variables according to 



u 



w 



du. (30) 



(31) 



1 + r' 1 + z' 

and integrating both side of (BH) with respect to d(f){x), taking into account (|29|), leads to 



"1— a— c 



JO 



[c + t{c — /x)] exp [—x{z — T)]dTd(p{x) . (32) 
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Taking the Laplace transform of the above identity we obtain for the Stieltjes transform 
of the measure d(j){x) the relationship 



s{p) 



d(f){x) 



^'(c+l,l-a;p) 



x + p \E'(c, — a;;p) + (c — yu)\E'(c + 1, 1 — ■ 



(33) 



which we can rewrite formaly, using the expression of the Tricomi function \E' in terms of 
generalized hypergeometric functions page 257], on the form 



s{p) =p"^(c+l,l-a;p)3Fi 



^c,c + a, ^ ^ ^ + 1 1 

c{c+a) ' ~ 



(34) 



where the principal branch of the 3F1 is considered as that one of the ^ function. The 
function ^(a, 6;p) having no zeros for | argp| < vr the denominator of (^) has no zeros in 
this region at least for fi < c. 

The CAL polynomials belong to the Laguerre-Hahn family of orthogonal polynomials and 
are of class zero ||22|. It is easily verified that the Stieltjes transform s{p) of the measure, 
calculated in (p3D, is a solution of the Riccati equation 

(35) 



ps'{p) = [fip - {c- fi){a + c- /i)] s^(p) + [p + a + 2{c - /i)] s{p) - 1. 



The absolutely continuous part of the measure c?0(x) can be computed using the Perron- 
Stieltjes inversion formula. Details of the method are in [M and we obtain 



0'(x) 



1 



r(c+ l)r(c + a + 1) |^(c, -a;xe''') + (c - /i)^'(c + 1, 1 - a;xe''^)|^' 
we can write formaly 



(36) 



0'(x) 



r(c+l,c + a + ll 



c, c + a, 

c(c+a) 



X 



(37) 



The CAL polynomials L'^{x]c,fi) satisfy the orthogonality relation, valid at least for 
fi<c, c>0,a + c> —1, 



L"(x; c, fi)L'^{x; c, fi)d(t){x) 



2.3 Fourth- order differential equation 



(c + a + l) 

(C+I)n 



{31 



The CAL polynomials verify a fourth-order differential equation p, ^ . One way to obtain 
this fourth-order differential equation is to start from their explicit form (^). The righthand 
side of ( |22| ) is a sum of two products of a iFi times a 2-^2- The iFi are solution of a second- 
order differential equation but for the 2-^2 a third-order one is expected. In fact the 2F2 
involved in (12^) are of the form 



y{x) 



2-^2 



X 



(39) 
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which can be shown with httle effort to be a solution of the second-order differential equation 

x[{b - e)x +e{d - e - l)]y"ix) - {(6 - e)x'^ + [e{2d - e - 2) + 6(1 - d)]x (40) 
-ed{d - e - l)K(x) - b[{b - e)x + {e + l){d - e - l)]yix) = 0. 

So we can use the Orr method to obtain the differential equation satisfied by the products 



involved in ( p2D p3| . Let us notice that the second product being obtained by the trans- 
formation T of the first one, the fourth-order differential equation will have to be invariant 
under this transformation. 

— c — a 



The function y 



X is solution of 



and the function z 



2-^2 



1 — a ' / 
xy" + (1 — a + x)y' + (1 + n + c)y 
— c, /i — c + 1 



(41) 



1 + a, — c 



:x \ of 



x[jjix +(c — /i)(c + a — ij)]z\x) — {fix^ + \{c — ij){c + a — ji) — a/ija; 
— {a + l)(c — /i)(c + a — /i)}z'(x) + c[{^x + (c — /i + l)(c + a — /i)]2;(a;) = 0. 

Changing the functions y and z to y = fv and z = with 

/ 



(42) 



Q — 1 X_ 

X 2 e2 . 
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[{c — ^){c + a — fi) + fix]2x 2^62, 



(43) 
(44) 



we obtain the normal form of the differential equations (|4l| ) and (^31) 

t;" + It; = 0, w" + Jw = 0. (45) 
The product u = vw is solution of the fourth-order differential equation (see page 146]) 



d_ 

dx 



u'" + 2{I + J)u' + {!' + J')u 
I -J 



-(/ - J)u 



(46) 



Finally we obtain the needed equation for the CAL polynomials setting y{x) = fgu. 

Details of this calculations are very difficult to write explicitly and were achieved with 
the help of the MAPLE computer algebra 0]. Although with this help the fourth-order 
differential equation for the CAL polynomials is not easy to find. We give it as a curiosity: 



C4y^^\x) 



with 



c^y^^^x) + C2y^^\x) + ciy^^\x) + coy{x) 



Ci=x\2Ax'^ + BX + 2C), 

C3 = 2x{3Ax^ + 2Bx + 5C), 

C2 = -2Ax^ + Dx^ + Ex^ + Fx + G, 

ci = -AAx"^ + Hx"^ - AAx^ + Ix + J, 

Co = n{n + l){2Ax'^ + Kx + 2L), 



0, 



(47) 

(48) 
(49) 
(50) 
(51) 
(52) 
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where 

A = fi^{l + 2n), 

B = /i(2(l + 4:n)f? - (4(1 + 2n){2c + a) - + 2c(3 + 4n)(c + a)), 
C = (c - ^)(c + « - /i)(2n/i^ - ((1 + 2?2)(2c + a) + + 2c(n + l)(c + a)), 
D = -/i(2(l + 4n)fi'^ - (2(1 + 2n)(8c + 4a + 1 + 2n) - l)/i + 2c(3 + 4n)(c + a)), 
E = _4^^4^((i + 4n)(4n + 12c + 6a + 3) + 3)/i^-2(l + 2n)((2c + a) 
X (4n + 12c + 6a + 2) - 2c(c + a) - l)/i^ + c(c + a) ((3 + 4n) 
X {An + 12c + 6a + 1) + 3)/i - 4c^(c + a)^(n + 1), 
F = 8n(2c + a + n + l)/ 

-((1 + 4n)((2c + a)(4n + 12c + 6a + 3) - 8c(c + a) - 1/2) + 10c + 5a + 5/2)^^ 
+ ((1 + 4n)((a^ + 6c(c + a)(2n + 8c + 4a + 3/2) - (2c + a)(12c(c + a) + 1/2)) 
+ {2a^ + 6c(c + a) - 1/4) (4c + 2a + 23/6) - 25/6a2 + 71/24)/!^ 
-2c(c + a)((l + 4n)((2c + a){2n + 4c + 2a + 3/2) + 2a^ - 1/4) 
+ (2c + a)(8c + 4a + 5/2) + 2a^ - 7/4)/i + 8c^(c + a)\n + l)(n + 2c + a), 
G = -2(a-2)(a + 2)(c-/i)(c+a-/i)(2(c-/i)(c + a-/i)n 

-(2c + a + l)/i + 2c(c + a)), 
H = 2/i(-2(5n + + ((1 + 2n)(n + 12c + 6a + 1/2) - 3/2)/i - 2c(c + a)(5n + 4)), 
/ = -Unn^ + 2{{n + l/5)(8n + 40c + 20a + 7/5) + 93/25)/i^ 

+2(-(l + 2n)((2c + a)(4n + 12c + 6a + 2) + 10c(c + a) - 1) - 6c - 3a)n^ 
+2c(c + a)((4/5 + n)(8n + 40c + 20a + 33/5) + 93/25)/i - 12c2(c + aY{n + 1), 
J = 4{fi-c){fi-c- a){3n{n + 2c + a + + (-1/4(1 + 2n){{2c + a){Qn + 8c 
+4a + 3) + 8c(c + a) + 8) - 9/2c - 9/4a)/i + 3c(c + a){n + l){n + 2c + a)), 
K = /i(2(4n - - (4(1 + 2n)(2c + a) - 3)fx + 2c(4n + 5)(c + a)), 
L = {c- fi){c + a- fi){2{n - l)fi^ - ((1 + 2n)(2c + a) + 6)/i + 2c{n + 2)(c + a)). 



Note the invariance of the differential equation (]47|) by the transformation T defined in (p!3|) 



2.4 Particular cases 

We now give the different results corresponding to limiting cases of special interest. 
2.4.1 Limit c = 

In this limit we obtain from (^) the co-recursive Laguerre polynomials. 



a 1 n\ \l-a,;U-a / \ 1 + a ' 
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The limit /i = in (^) leads back to the classical Laguerre polynomials. An explicit form is 

(a + 1) " (-n)fc 



^n(^;0,/i) 



fc=0-'(l + «)fc 

fi{k — n) 



1 + 



l + k){l + a + kj^ \ k + a + 2,k + 2 



1 + k — n,l + a,l ^ 



x\ (54) 



and the corresponding absolutely continuous part of the measure is given for fi < 0, a > — 1, 

by 

-2 



T(l + a] 



1 — 1 — a;xe 



(55) 



where the limit /i = is also straightforward. 

It is easy to see that the differential equation (^) satisfied by the co-recursive Laguerre 
polynomials can be factorized in the limit c = to obtain the fourth-order factorized (2+2) 
differential equation 

xAix)!)^ + {(2 + a - x)A{x) - xB{x)}D + {n - l)A{x) + {x - a - l)B{x) + C{x) 

X [xD^ + {x + l-a)D + n+l\ L^(x; 0, fx) = 0, (56) 

where D = d/dx and 

A{x) = 3x + 2(x — a + fi){2n{x — a + /i) + x — a — 1}, (57) 
B{x) = 1 + 2x - 2a + 2(1 + 4n)(a; - a + fi), (58) 
C{x) = (1 + 2a; - 2a){l + a-x- 2n{x - a + fi)} + 3x{l + An). (59) 

The comparison with the differential equation given in |2^, Eq. 34-35] requires some atten- 
tion because of a few misprints. 

2.4.2 Limit c = —a 

In this limit we obtain a special class of CAL polynomials corresponding to the associated 
Laguerre polynomials for which 



L^{x;-a) 



n\ 



[I - a)n 



L-'^ix). 



(60) 



We can write the L"{x] —a,fi) 

rai \ e"^ f (-a - /i)(l - a)„ + a + 1 \ ( -n 

[l — a)n—a [ n\ \^i + a,fi + a J — 

„ f—n + a \ I 



a 



X 



(61) 



which gives ( |60D in the limit /i = 0. Except for the global factor jiz^, ( pTf ) is obtained from 
(p3D changing a to —a. The corresponding measure and differential equation are obtained 
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in the same way from p.4.1| . An explicit form is 



L^{x; -a,fi) = }^ — 



1 + 



fi{k — n) 



p p + A;-n, l-a,l_^ 



;i + A;)(l -a + A;)' '\k~a + 2,k + 2 



X 



(62) 



2.4.3 Limit fi = 



In this hmit we obtain the associated Laguerre polynomials studied in and |T2| 
with the measure 



^ n , T-N(c + a)n+i (1-c-a \ f-n-c 
— (1 + T) iFi 1 ;a; iFi . ;x 



[X 



An explicit form is 



L°(x;c) 



r(l + c, 1 + c + a) 
c + a + 1)„ " 



^'(c, 1 - a;xe~''') 



(63) 



(64) 



E 



n\ ^ (c+ l)fc(c + a + l)k 

k — n,c,c + a 



X3F2 



.c + k + l,c + a + k + l 



a\x' 



(65) 



The limit c = leads back to the Laguerre polynomial case. The coefficients of the differential 
equation (^) satisfied by the associated Laguerre polynomials are now very simple 

C4 = x^, C3 = 5x, C2 = —x{x — 2F) — + 4, ci = 3(F — x), cq = n{n + 2), (66) 

with F = n + 2c + a. This differential equation was first given by Hahn Eq. 22]. See 
for the special factorizable case c = 1 and ^ when c is an integer. 



2.4.4 Limit fi = c 



In this limit we obtain the so called zero related Laguerre polynomials studied in |T2|. Note 
the symmetry T of the monic polynomials is now broken. 



^"(^,c) = e- 



and the measure 



[X 



(c + a + 1) 
(c+l)„ 

a{a + 1) \2 + a \ I- a 



\E'(c, —a; xe 



r(l + c, 1 + c + a 
Again the limit c = leads back to the Laguerre polynomial case. 



(67) 



(68) 
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The explicit form ( [25| ) simplify in 



L:(.;c.,)=(i±^t 



k — n,c,c + a + l 



^- ^Q{c+l)k[c + a + l)k 



The coefficients of the differential equation satisfied by the zero related Laguerre poly- 
nomials are 



C4 = x^{2{2n + l)x + D), C3 = 2a;(3(2n + l)x + 2D), 

C2 = -2(2n + l)x^ + (8n(F + 1) + 8c + 

-(4(a^ - l)n - 2o? - D(4c+ 1) - l)x - 1/4L'(L)2 - 9), 
ci = 8(2?2 + \)x^ - 4(2n(F + 1) + 2c - L))a; - 2D{n[D + 3) + 6c + 2D), 
Co = n[n + l)(2(2n + \)x + 3D), 



with 



F = n + 2c + Q!, D = l + 2a, 
and are no longer invariant under the transformation T. 



(70) 
(71) 

(72) 
(73) 



(74) 



2.4.5 Limit = c + a 



This is a new simple case of CAL polynomials lacking in |T2]. 



£^(x,c) = e-^|,F,^-;;.j,F,^-7_^;";.j 

I (c + a)(c + a + l)„ /i-c-a \ /-^-c 
a(l-a)(c+l)„ \2-« / + a 



and the measure is obtained using (10) page 258] 

^'(c+ 1,2 - a\xe 



5 I f 5 



r(l + c,l + c + a) 



(75) 



(76) 



In this case the limit c = does not lead back to the Laguerre polynomial case but to the 
co-recursive Laguerre one with fi = a. 
The explicit form is 



k — n,c + l,c + a 



(c + a + l)n " , 

n\ toi^ + l)fc(c + a + l)r'[c + k + l,c + a + k + 



1 ; l) x'. (77) 



The coefficients of the differential equation (|^) satisfied by the polynomials £"(a;, c) are 
obtained from ([fO|-[7^) changing only D = 1 + 2q; by D = 1 — 2a. 
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3 The case of Jacobi polynomials 

We now present some results on the CAJ polynomials. The recurrence relation of the asso- 
ciated Jacobi polynomials P^'^{x; c) is ^1 



{2n + 2c + a + p + 1) [{2n + 2c + a + p + 2){2n + 2c + a + p)x + - p^]pn 

= 2{n + c + l){n + c + a + p + l){2n + 2c + a + P)pn+i (78) 
+2{n + c + a){n + c + p){2n + 2c + a + p + 2)pn-i- 

We can note the invariance of the recurrence relation (|79D under the transformation T' 
defined by 

r'(c,a,/3) = (c + a + /3, (79) 

All polynomials satisfying (^), for which the initial conditions are symmetric in a and P 
and invariant under T' have the property 



p;"'-^(x,c + a + /3)=K'^(x,c). 
As in [^1 we use the more convenient shifted polynomials defined as 

R''/{x;c) = P,^''^{2x -l;c). 
Due also to the properties of the recurrence relation (79) we have 

R?-'f(x- c) = (-l)"/2f'°ri-x;c). 



^0) 



^2) 



3.1 Explicit representation for the CAJ polynomials 

A solution of the recurrence relation satisfied by the R"'^{x; c) in terms of the hypergeometric 
function is ET], page 280] 



Ur 



c + a + l)n „ (-n-c,n + c + a + P + l . ^ 
2Fi\ , , ; 1 - a; 



(c+l)„ ' '\ l + « 

and an other linearly independent solution is given by 

^, (c + /? + !)„ f-n-c-a- p,n + c+l . \ 

Vn = T Un = 7 ■ --5— il~2-Fl . ] 1 - X ] . 

{c + a + p + l)n \ J 



^3) 



M) 



The functions m„ and x ^(1 — x) ""Vn are two independent solutions of the second-order 
differential equation 

x{l - x)y"{x) + [l + P-{a + p + 2)x]y'{x) + {n + c){n + c + a + p + l)y{x) = 0. (85) 
The associated Jacobi polynomials are defined by (|79|) and the initial condition 

P-f(x;c)=0, Po^'^{x;c) = l. (86) 
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This gives for P"'^(x; c 



_^ _ {2c + a + (3 + l){2c + a + (3 + 2) 



2{c+l){c + a + (5 + 1) 

«2 _ ^2 



X + 



{2c + a + (3){2c + a + (3 + 2) 



17) 



The CAJ polynomials P^'^ {x; c, fi) satisfy the recurrence relation ([79|) with a shift on 
the first monic polynomial. This corresponds to the initial condition on the shifted CAJ 
polynomials c, /i) 

it_i (x,c,/ij - 2(c + q;)(c + /3) ^ ° l^;, c,/ij-l. 

Ifc + a;^0orc + /5^0 this initial condition in (^) leads nevertheless to a shift yU on the 
value of X in '^(x; c). 
As in section writing 

P^'^(x;c,/i) = AM, + Pt;„ (89) 



and using (|88|) we obtain 

= 1 [D^;o - t^-i] and P = [Duq - , (90) 

where A is easily calculated using the fact that m„ and x~^{l — x)~°'Vn are two independent 
solutions of (|85D , 

a(2c + a + /5) 
(c + a)(c + /?) 

The condition A 7^ leads to a 7^ and 2c + a + (3 7^ 0. We can note the invariance of D 
under T' and that B = T'A. 

Grouping the two 2P1 involved in the expression (^O]) of A gives 



c + a + P-D{c + P) /-c-a-(3,c,F + l \ 

ai2c + a + (3) [ l-a,P ' ^ " ^ J ' ^^2) 



(93) 



with 

c[P(c + /3)-c-a-/3] 

P(c + ;3)+C 

and the CAJ polynomials could be writed 

fir (:^; = (1 + TO ° + ° + ^-^';;^\ c + 

a(2c + a + /?) (c+l)„ 

„ /-c - a - /5,c, P + 1 ^ \ _ f-n-c,n + c + a + (3 + 1 ^ \ 



We will use this expression of the CAJ polynomials in p.4| to obtain a fourth-order differential 
equation of them. 
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Transforming the 2-^1 (1 — x) in ( ]89| ) by Eq. 1, page 108] one obtains with httle algebra 

4)Mc + a) 



P{2c + a + p) 



{c + a + l)n 
'c + a + (3 + l)n 



—n 



a — /?,n + c+ l 



X 



(95) 



-c,c + a + [3 + 1 
1 + 13 



This formula generalizes the one of Eq. 28] to the case of the CAJ polynomials. As the 
explicit form of the CAL polynomials (^) the representation (0) and (|95|) are valid only 
for a 7^ 0, ±1, ±2 . . . and /3 7^ 0, ±1, ±2 . . . but can be extended by limiting processes. 

We obtain an explicit formula following the same way as in We first use Eq. 14, 
page 87] for each product of 2F1 in (|95D to obtain four series involving gamma functions and 
a 4F3. For two of them we use Eq. 1, page 56]. The next step is to use for each 4F3 
twice 0, Eq. 3, page 62]. After numerous cancellations only two series of 4F3 remains we 
can group to obtain the following explicit form 

Da,/3/ . ^ _ r_i Nn (2c + a + /3 + l)n(/3 + c+l)n A (-n)fc(n + 2c + a + /3 + l)fc 



n\{c + a + 13 + 1) 



k=0 



(c+l)fc(c + /?+l)fc 



X5F4 



where 



G 



^k-n,n + k + 2c + a + p + l,c,c + P,G + l 
^ c + k + l,c + p + k + l,2c + a + P + l,G 

2c{c + f3){2c + a + (3) 



2c(c + (3)+ /i(2c + a + I3){2c + a + (3 + 1)' 



1 x^(96) 



(97) 



3.2 Generating function 

One can obtain a generating function of the CAJ polynomials following the same strategy 
pl[] for the associated one. Let Q{x, w) be a generating function of i?^'^(x; c, /i) 



as m 



N ^{c+l)nic + a + P + l)n , 

nli2c + a + P + 2U " 
Starting from the form ( p4|) for the i?"'^(x; c, /x) it follows 



(9^ 



rf ^ n . ^^^ c + a + P-Djc + P) ^ 

^(-' -) = (1 + r ) a(2c + a + /3) + "^^"^^ 

^ (c + a + l)„(c + a + /3+l)„ „ /- 
^ „i^o -!(2c + a + /5 + 2). " [ 

and using |^ Th. 4] we obtain 

rf ^ n , ^;, c + a + P-Dic + (3) 



c — a — /3, c, F + 1 
1^ l-a,F ' 

-n — c,n + c + a + (3 + 1 
l + a 



X 



1 - X , (99) 
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where 



X (c + a) 



W{Z2+1)_ 



(-c~a- p,c,F + 1 \ 



(100) 



(-c,c + a + (3 + 1 I- ZA /c + a + l,c + a + (3 + 1 2 
M l + a ' 2 T M 2c + « + /5 + 2 'I + Z2 



^1 



1 



1 + wY — Awx 



l + J(l + w) 



Awx 



w w 
which generahze the aheady exotic generating function Eq. 75]. 



(101) 



3.3 Spectral measure 

The Stieltjes transform of the measure of the shifted associated Jacobi polynomials R"'^{x; c) 
is iH, Eq. 63-64] 



. . 1 P /c+l,c + ^ + l 1 , „ 
s{p) = -2i^i I o„ , „ , ^ , o ; - 1 2i^i 

p 



2c + a + /3 + 2 'p 



^ c,c + f3 r 

.2c + a + (3' p 



(102) 



The CAJ polynomials i?"'^(x; c, /i) satisfy the same recurrence relations as the -R"'^(x; c) 
with a shift fi on the first monic polynomials 



i?i'^(x; c, /i) — i?i''^(x; c 



(2c + a + /3 + l)(2c + a + /3 + 2) 
2(c+l)(c + a + /5 + l) 



-/i. 



(103) 



Using continued J- fractions |T^, ||, ^ whose denominators are i?°'^(x; c, /i) and i?°'^(x; c) 
we can derive for the Stieltjes transform of the measure of the CAJ polynomials 



s(p; /i) = s{p) { 1 + f^s{p) 



^c+l,c + p + l_V 
y2c + a + (3 + 2' p 



c,c + (3 .n /i /c+l,c + /5 + l.l 
' \2c + a + (3' py 2^ '\2c + a + (3 + 2' p 



(104) 



that we can also write using contiguous relations 

^c + l,c + P + 1 1\ 



sip; At) = 2F1 



2c + a + (3 + 2' pj 

c + a 2c + a + p+l 



+ 



2c + a + (3 

c + (3 



+ 



2c 

2c + a + (3 + 1 



2c + a + (3 ' 2c 
Grouping the 2-^1 in ( p.05| ) gives the compact formula 

^c+l,c + (3 + 1 1 



/i) 2F1 



c,c + (3 . 1 \ 
2c + a + (3 + 1' p) 

' c,c + p + l _V 
2c + a + (3 + 1' p 



(105) 



s{p; /i) = 2F1 



2c + a + (3 + 2' p 



3F2 



c,c + (3,G + i r 

2c + a + (3 + l,G' p 



(106) 
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where G is given by (0). A sufficient condition for the positivity of the denominator in 
( |106D on (l,oo) is 



c > 0, 



c> 



a > — 1, 



fx 



> 



2c{c + f3) 



{2c + a + I3){2c + a + (3 + 1)' 



(107) 



but other conditions are possible. 

To obtain the absolutely continuous part of the spectral measure we need to evaluate 
s~^{p]fi) — s^{p;fi) where are the values of s above and below the cut [0,1]. Using the 
analytic continuation 0, Eq. 2, page 108] for each 2-^1 in ( |104D we find for the spectral 
measure of the CAJ polynomials 



= (1 
= (1 - 



c,c + P _e^\ n c+l,c + P + l_e^ 
' M2c + a + /? ' X i 2x ' 'Uc + « + /5 + 2 ' x 



3-^2 



c,c + p,G + l e* 



:i08) 



^2c + a + p + l,G' X 
valid at least under the conditions (|107|) . 

3.4 Fourth- order differential equation 

The method used to obtain the differential equation satisfied by the c, /i) is the same 



as in ^73| . In (|9^) the hypergeometric function 2-^1 is solution of the equation (|85D and the 
3F2 is of the form ^F^, ^■'^■'^ 
equation 



qx-o 1 ' ', \x\ which is also solution of the second-order differential 



x{x — 1) [(a — e)(6 — e)x + e{d — e — 1)] y"{x) 

+ {(a -e){h-e){a + h + l)x^ + [e(a + 6 + l)(2rf - e - 2) - d(a6 + e^) + afe] a; (109) 
+ de{e -d + 1)} i/'(x) + ab [(a - e)(6 - e)x + (e + l)(rf - e - 1)] y{x) = 0. 

We don't write here the fourth-order differential equation hardly obtained by symbolic 
MAPLE computation. The coefficients are at most of degree eight in x and it would take 
several pages to write them. We give the results only in the following limiting cases. 



3.5 Particular cases 
3.5.1 Laguerre case limit 

The limit giving the CAL polynomial case is obtained by the replacement 



X 



1 _ 2x 

^ /3 
,2^ 



P 



00, 



(110) 



in P"'^(x;c, fi). The representation (|9^) is the more suitable to obtain the form of the 
CAL polynomials (53) using the Kummer's transformation (pSf) for one of the confluent 
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hypergeometric functions and his generahzation 



/a,e+ 1 



X 



'c — a — 1 



e(c— a— 1) 



e—a 
e(c— a— 1) 



+ 1 



-X 



(111) 



for one of the 2-^2- Note that ( |1 1 1| ) gives ( |2BD in the hmit e ^ 00. 



3.5.2 Limit c = 

In this hmit we obtain the co-recursive Jacobi polynomials. An explicit form is 



K^{x-fi) 



k=0 



(112) 



/ fx{k-n){n + k + a + f3 + l) fk - n + l,n + k + a + p + 2, P + 1,1 \\ 
^\ 2{k + 1){(3 + k + 1) k + 2,(3 + k + 2,a + f3 + 2 ' ^ j J ' 

and the spectral measure is given by 
(f)'{x) = (1 - x^xf^ 



1,P + 1 e' 



1 ^ ii F 

' la + /3 + 2' X 



-2 



(113) 



The limit /i = leads back to the Jacobi polynomials. 

The fourth-order differential equation satisfied by the co-recursive Laguerre polynomials 
can be factorized in the limit c = to obtain as in ||2^ the factorized (2+2) differential 
equation 

= [(1 - x^)A{x)D'^ + {(/?-« - (« + /? + 4)x)A{x) - (1 - x'^)B{x)]d 
+ [n{n + a + p + l)-{a + P + 2)} A(x) + {/?-«-(« + /? + 2)x} B{x) + C{x)] (114) 
X "(1 - x^)D2 + {(a + P-2)x + a- P}D + n{n + a + p + 1) + a + p] K'^(x; /i), 

where 

A(x) = 2{a + PY{2n + l){n + a + p + l/2)x'^ + 2{a + /?)(4n(n + a + P + 1) 
X (-/i(l + a + /3) + a-/3) + (l + a + /5)(-/i(a + /5 + 2) + 2a - 2p))x 
+4n{n + a + P + l)(-/i(l + a + /5) + a - /3)^ - (a + /3)(-2/i(l + a + /3) 
X (/5 - a) - 2(/3 - af - 3a - 3/5) 

5(x) = -(a + /5)((a + /3)(8n(n + a + p + l) + 3a + 3P)x + 8n{n + a + p + 1) 

x(-/i(l + a + p) + a - P) - 2{a + p + 2){1 + a + P)fi + {a - p){3a + 3P + 4)) 

^(x) = -(a + /5)((a + /3)(a + /5 + 2)(2n(n + a + p + l) + a + P- l)x^ 

+2(n(n + a + /3 + l)(-/i(a + P + l){a + P - A) + 2{a - P){a + P - 2)) 
+ (a + /?)(a - /5)(a + - l))x - 2n(n + a + /5 + l)(-/i(a + /5 + 1)(/? - a) 
-(a - py + 6a + 6/3) + (a + /3)((a - /?)2 - 3a - 3/3 - 6)) 
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3.5.3 Limit c = —a — f3 

Due to the T' invariance of (|79|) we obtain in tliis limit the special case of CAJ polynomials 
for which 

R^'^ix; -a - = (115) 
All the results are obtained from |3.5.2| by changing a to —a and j3 to —f3. 



3.5.4 Limit c = -/? 

The explicit form ( p6D simplifies in the same way as in the case c = 0. One obtain 

KAx; -ft A) - (-1)" '7 i t ! ° :f ^ ^'^ x^ (116) 

(a + l)n ^0 fc!(l-/?)fc 
[ /x(A;-n)(n + A; + a-/? + l) /fc - n + 1, n + + a - /5 + 2, 1 - /3, 1 \\ 
""{^ 2{k + - [3 + k) ' '1^ A; + 2,2-/5 + fc,a-/3 + 2 ' ^ j / " 



Comparing this form with the explicit form of the co-recursive Jacobi polynomials ( |112D one 



see 

R-J{x; -P, /i) = .^'^""f/^^; /.), (117) 

of course the spectral measure and the fourth-order differential equation are obtained from 
( p2D and ( p4|) changing /3 to 



3.5.5 Limit c = —a 



This case is the T' transform of the preceding case. All the results are obtained from |3.5.4 
by changing a to —a and P to — /?. 

3.5.6 Limit ^ = 



In this limit we obtain the associated Jacobi polynomials studied in [^. The form [^T], Eq. 
28] is obtain directly using (|95|) but a form slightly different is 



?a,/3r™. „N _ n , -^/A {c + a){c + a + I3){c + a + 1), 



i?r(^;c) = (i + T')- 



a(2c + a-F/3) (c-M) 



„/l-c-a-/9,c^ \ _ /-n - c,n + c + a + /? + 1 . \ 
x,FJ ;1-xUfJ ;l-x , (118) 



The explicit form ||HT], Eq. 19] is easily obtain starting from ( PBf ) with G = 2c + a -|- /3, the 
5F4 reducing to a 4F3. Obviously the limit c = leads back to the Jacobi polynomials. 

The coefficients of the differential equation ( ^7|) satisfied by the associated Jacobi poly- 
nomials are 

C4 = x^{x-lf, C3 = 5x{x-l){2x-l), (119) 

C2 = (24 - (ra + 1)^ - a) 2;(a; - 1) - 5x - /^^ + 4, (120) 

Cl = -3/2 ((3A + (n + 3)(n - l))(2x - 1) + 5) , (121) 
co = n(n + 2)A, (122) 
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with 



A = (C + n + l){C + n-l), B = {a - P){a + P), C = 2c + a + p. 



(123) 



This result was also first given by Hahn Eq. 20]. Note the T' invariance of A, B and C 
leading to the invariance of the q more obvious than in 1^1], Eq. 47-48]. 

3.5.7 Limit ^ = (gc+zj+^/S+ij+a+i) 

In this limit the symmetry T' is broken. We obtain the zero-related Jacobi polynomials 



studied in |jT5[ . An explicit form is 

(2c + a + /3 +!)„(/? + c + 1) 



K'^(x;c) = (-iy 



y 

n\{c + a + P + fr'o {c+l)k{c + p + l)k 

k — n,n + k + 2c + a + f3 + l,c, c + f3+l \ ^ 
c + k + l,c + p + k + l,2c + a + P + l '^)^- 



-n)kin + 2c + a + P + 1) 



(124) 



The limit c = leads back to the Jacobi polynomials and the limit defined in ( |110| ) gives the 
zero-related Laguerre polynomials ( |2.4.4|) , using Kummer's transformations. The spectral 
measure is 

^ c,c + /5 + l _e^^ 
^2c + a + /3 + 1 ' X ^ 

The coefficients of the differential equation (^71) satisfied by the polynomials 7^°'^(x;c) 



0'(x) 



'1 



(125) 



are 



C4 



C2 



Cl 



x2(x- l)2(Ax + L)), (126) 

C3 = x(a; - 1) {sAx^ - 3{A - 3D)x - AD) , (127) 

-1/2A{A + 2^2 - 29)x^ + (l/2A{A + 2^^ - 2fi - 23) - D{C^ - 19)) (128) 
-1/4 (A(D + 1){D - 3) - 2D(2(:72 _ 2fi + D - 35)) x - 1/AD{D^ - 9), 

-A{A + 2C^ - 5)x^ + 1/4 (A(A + 2^2 - 2B - 5D - 5) - 3D{AC^ - 11)) a; (129) 
+ 1/4D ((D + 3)A + QC^ -QB + 3D- 15) , 

2n{n + 1){C + n){C + n + l){Ax + 3D), (130) 



Co 



where B and C are defined in (11231) and 



A = {2n+l){l + 2C + 2n), 

2{c+l3){c+a+l3) 



D 



2/3. 



(131) 



3.5.8 Limit - (^^^+2c)(/3+a+2c+i) 

This case is the T' transform of the case |3.5.7| . The explicit form is 

{-n)k{n + 2c + a + P + l)k 
k=o {c + a + (3 + l)k{c + a + l)k 



X4F3 



n!(c+l)„ 

k — n,n + k + 2c + a + (3 + l,c + a + (3,c + a + 1 
c + a + (3 + k + l,c + a + k + l,2c + a + (3 + 1 



■l\x^. (132) 
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The limit ( [110| ) leads back to the Laguerre case p.4.5| and the limit c = to the co-recursive 
Jacobi polynomials with /i = ■ The spectral measure is 



b'{x) = (1 - 



' c+l,c + /? .e^ 
,2c + a + /5 + l' a; 



(133) 



The coefficients of the differential equation (^7|) satisfied by the polynomials lK"''^(a;;c) 



are obtained from (|T26|-|T3T|) changing only D = 1 + 2/5 by D = 1 — 2/3. 



3.5.9 Limit /x 



2c(c+/3) 



(2c+a+/3){2c+a+/3+l) 



The symmetry T' is also broken. We obtain a new simple case of CAJ polynomials. An 
explicit form is 



= { ^^„ (2c + a + /3 + l)n(/3 + c+l)„ ^ (-n)fc(n + 2c + a + /3 + 1)^ 



n!(c + a + /? + 1)„ 



fc=0 



(c+l)fc(c + /5 + l)^ 



/fc-n, n + A; + 2c + a + /5+l,c, c + /? \ ^ 
Ic + A; + 1, c + /5 + + 1, 2c + a + /5 + 1 ' i ^ ■ 



(134) 



and the spectral measure 



0'(x) 



:i - x)"x^+2c 



c,c + (3 
2c + a + /3 + 1 



(135) 



The coefficients of the differential equation (^7|) satisfied by the polynomials 7^°'^(a;;c) 



are 



C4 = x2(x - 1)2 (A(x - 1) - D) , 
C3 = x{x - 1) (sAx^ - {ISA + 9D)x + 5A + 5D 

C2 = -1/2A{A + 2C^ - 29)x^ + (^A{A + 2C^-B- 32) + D{C^ - 19) 

-1/2 (A(A + 2C2 + 2/3^ - 25 - 43) + D(2C2 -2B-D- 41)) a; 

+(/52-4)(A + D), 
ci = -2A{A + 2^2 - 5)^2 

+ 1/2 (A(7A + 14^2 -2B + 5D- 35) + D(3C2 - 33)) x 

-A{3/2A + 3C^ - B- 2a^ - 7) - ?,D{C^ -B- 1/2D - 3), 
Co = n(n + 1)(C + n)(C + n + l){A{x - 1) - 3^), 



(136) 
(137) 

(138) 



(139) 



(140) 



where B and C are defined in (D-231) and 



A= {2n + l){l + 2C + 2n), D = l + 2a. 



(141) 
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C5 c in TlmH//- 2(c+a)(c+a+/3) 
Ci.b.iU L,imit /i - - (2c+a+/3)(2c+a+/3+l) 

This case is the T' transform of the case |3.5.9 . The exphcit form is 



1 = r-1 ^n (2c+a + /^+l)n(ft + C+l)n " (-n)fc (n + 2c + ft + + l)fc 

"^"''''^ ^ ^ n\{c+l)n t^^{c + a + (3 + lUc + a + l)k 

(142) 



X4F3 



fc=0 

''k — n,n + k + 2c + a + P + l,c + a + (3,c + a 
.c + a + p + k + l,c + a + k + l,2c + a + p + 



and the spectral measure 



= (1 - ^^"-2^/3+2c+2 



^c+l,c + /3 + l _ 
i2c + a + /5+ 1 



X 



(143) 



The coefficients of the differential equation ( ^T] ) satisfied by the polynomials 9^"'^(x;c) 
are obtained from (|136Hl41| ) by changing only D = 1 + 2ahy D = 1 — 2a. 

3.6 Conclusion 

We end by brief remarks. In this article we have studied properties of the co-recursive 
associated Laguerre and Jacobi polynomials which are of interest in the resolution of some 
birth and death processes with and without absorption. For few values of the co-recursivity 
parameter we obtain polynomial families for which the results are of the same complexity 
as the corresponding associated polynomials. For the CAL polynomials we find the two 
expected cases corresponding to = /^o (zero-related polynomials) and the new dual case 
/i = A_i For the CAJ polynomials, due to the properties ( ^OD and (|^) we have two 

more cases corresponding to fi = —T"fio and fi = — T"A_i where the transformation T" is 
defined by 

r"{c,a,p) = {c + a + p,-p,-a). (144) 

In some cases the fourth-order differential equations satisfied by the polynomials studied 
above are factorizable (co-recursive and associated of order one) but we don't find factoriza- 
tion either for the co-recursive associated polynomials or for the associated one. Of course 
this is not a proof that the conjectures on this factorizability made in are wrong. 
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